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Abstract 

We apply a number of schemes which variationally improve perturbation 
theory for the Kardar-Parisi-Zhang equation in order to extract estimates 
for the dynamic exponent z. The results for the various schemes show the 
same broad features, giving closer agreement with numerical simulations in 
low dimensions than self-consistent methods. They do, however, continue to 
predict that z = 2 in some critical dimension d c in disagreement with the 
findings of simulations. 
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The Kardar-Parisi-Zhang (KPZ) equation |I[ is perhaps the simplest nonlinear stochas- 
tic equation of the diffusion type. Nevertheless, the large-distance scaling behavior and 
anomalous dimensions have not yet been understood on the basis of systematic calcula- 
tional schemes, such as the renormalization group. To date the most reliable estimates for 
the anomalous dimensions, obtained by direct analytic means from the KPZ equation, have 
come from self-consistent or pseudo- variational procedures P-fl. These approaches are, in 
general, rather ad hoc and it would be useful to have a method of improving estimates 
and perhaps gaining some idea of their accuracy. The relative merits of procedures such as 
these, genuine variational procedures (by which we mean those with an associated bound) 
and "improvement" methods such as "the principle of minimal sensitivity" (PMS) ||, have 
recently been compared for stochastic processes described by simple Langevin equations (or 
equivalently Fokker-Planck equations) ||. In this Letter we extend these considerations to 
the KPZ equation — which can be formulated as a functional Fokker-Planck equation. It 
turns out that, while genuine variational techniques are not straightforward to use in this 
case, the application of the PMS yields improved values for the dynamic exponent z which 
agree better with the results of numerical simulation of models believed to be in the KPZ 
universality class. 

The KPZ equation for surface growth in (d + 1) dimensions with random deposition is 

h(x,t) = u V 2 h + g(Vh) 2 +v(x,t), (1) 

where the single-valued function h(5c, t) represents the interface and x is a rf-dimensional 
vector. The subscript "0" on surface tension Uq and noise strength Dq (below) is used to 
distinguish these bare parameters from the effective (renormalized) ones to be introduced 
later. The noise rj(x.,t) is Gaussian-distributed with zero mean (r/(x, t)) = and delta- 
function correlations 

(rifrtWmj)) = 2 J D 5(x-x / )(5(^-^ , )• (2) 
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Equivalently the noise may be specified in terms of the probability density functional: 



V[rj\ ~ expj-^A))- 1 Jdtjd d x[r)(2,t)f 



(3) 



The KPZ equation ([1]) can then be taken as a transformation of variables {r/(x, t) — > fo(x, £)} 
0] which leads to the following probability distribution of the h field: 



V[h] 



Det 



dh(5c, t) 



drj(-x!,t') 



exp 



4-Dn 



dt / cTa; 



/i - z/ V 2 /i - ^(V/i) 2 



(4) 



The Jacobian above can be shown to be h- independent; it plays no role in what follows and 
will be dropped. 

Introducing the response field /i(x, t) || yields the distribution of the h and h fields as: 
V[h, h] ~ exp{So[h, h]}, where 



S [h,h] = JdtJ d d x{D ~h 2 {x,t) - [h-v V 2 h- g{Vhf \ h(*,t)}. 

Note that the /i's are actually imaginary. Fourier transforming the fields gives 

~ f du f d d k 

S [h, h\ 



(5) 



2tt J (2n) d 
du r d d k r duj 



D h(k, uj)h{— k, — u) — [f^k — iooj h(k, oo)h(— k, — to] 

<J ' tk ' -[k ■ k']/i(k, u/)/i(k', w')/i(-k - k', -w - w 7 )- (6) 



2tt 7 (2vr) a! 7 2tt 7 (2vr) dl 

The problem is completely solvable in the absence of the nonlinear term (g = 0) |J, but 
one cannot access the strong-coupling fixed point through perturbation and the renormal- 
ization group. With these approaches ruled out, we turn to approximation techniques such 
as variational methods and related procedures. But just what kind of schemes are available? 

Before answering this question, let us make the first step toward a variational scheme by 
adding and subtracting an effective action functional S e s [h, h] 



V[h,h] ~ exp {S eS + A (<S - <S eS )} 



(7) 



where A serves as a counting device and will eventually be set to one. The most satisfactory 
approach would be a genuine variational scheme. Many of these are based on some version 
of the inequality e x > 1 + x which can be extended to 



2N-1 X n 



" > (8) 



n=0 



It is tempting to expand exp{A(iS>o — S c s)} and invoke this inequality; however, it does not 
hold as the h's are imaginary and the above inequality applies to real x only. 

Since the existence of an imaginary part to Sq\H, h] appears to be the stumbling block, 
perhaps it would be better to tackle the problem prior to introducing the response field. 
The inequality would hold, but a new problem arises. While one can proceed with a pertur- 
bation expansion around the bare quadratic terms in (f|), introducing variational quadratic 
terms to serve as a basis for perturbation rearranges the series, and one loses the property 
that divergences associated with the disconnected diagrams cancel order-by-order. These 
divergences render the inequality useless. 

There are other genuine variational schemes based on the differential form of the Fokker- 
Planck equation || and applicable to problems with a finite number of degrees of freedom. 
But they too are plagued with divergence problems when applied to a field theory. On 
the other hand, techniques are available which involve similar computational steps to the 
variational one, although the philosophy is rather different. In particular, they do not involve 
bounds such as that given by (H). These techniques can be formulated by choosing a form 
for the effective action S e g which is (i) quadratic (for calculational convenience) and (ii) 
has the same structure as quadratic terms in the bare action So, but with i>oq 2 — > z/ q and 
D — > Dq, that is, with the bare surface tension and bare noise strength replaced by their 
effective, or renormalized, counterparts: 

S e g[h,h] = J — J ^-p[Z) k /i(k,o;)/i(-k, -u) - (z/ k - iu) h(k, u)h(-k, -w)J. (9) 

The simplest scheme is probably the one pioneered by Edwards where one chooses 

the functions z/ k and -D k to be the "exact" (but as yet unknown) surface tension and noise 
strength respectively. When V[h,h] is expanded in A, the zeroth-order term will give the 
correct result (by construction) and hence the rest of the perturbation expansion (terms 
in A, A 2 , . . .) will have to be identically zero. This together with a scaling ansatz for the 



response and correlation functions, give values for the exponent z ||, which we will later 
compare to our findings. This scheme has been called "fastest apparent convergence" or 
FAC. 

The approach based on the PMS is calculationally similar to those already described 
- it involves the same diagrams, since one is expanding the same quantity about the same 
zeroth-order form (|9|), but now the rationale is different. The response or correlation function 
being calculated should not depend on i/ q or _D q , since these were introduced artificially 
and are not part of the "real" problem; however, any truncated expansion does depend on 
them. One then imposes a stationarity condition on the expansion to determine the result 
which displays the least dependence on them. Thus one asks that the expansion mimic as 
best as possible one feature of the true solution — its insensitivity to these variational terms 
— the hope being that other features will be mimicked as well. 

As we have stressed, all the schemes discussed so far are based on a perturbative calcu- 
lation in the parameter A. The structure of the perturbation expansion and the diagrams is 
essentially that given by Forster, Nelson and Stephen |TTJ], except that these authors were 
expanding in g and we are expanding in A. Thus we have quadratic, as well as cubic, terms 
in our "interaction" term and hence have terms odd in A in addition to the more familiar 
0(A 2 ) terms of Ref. |l 



A perturbative calculation of (\h(q, uj)\ 2 ) to 0(A 2 ) yields: 
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k • q 
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(k + q) 


' D 



(z/ q - iuo) (z/ q + iujf J (2n) d u k (u k + z/ k+q + iu) 
where the internal frequency integral has been carried out. 



+ c.c. 



(10) 



There are various ways to proceed at this stage. We want the calculation to be sim- 
ple without sacrificing any physics, so let us focus on the large-distance scaling behav- 
ior. Response and correlation functions calculated from (^) would be invariant under 
the scale transformation {x — > 6x, t — ► b z t, h — > b x h} if the parameters scaled like 



LZ-2 



V. 



D — > b d 2x+z D,g — > 6 2+x 2 g}. We build these scalings into our effective 



surface tension and noise strength by choosing 

Aq z and _D q 



Vex 



Bq 



-d-2 X +z 



(11) 



We also impose the exponent relation z + x — 2 due to the Galilean invariance fll|jl2|j . In 
what follows we will retain only the leading scaling behavior, thus the bare terms u and D 
will no longer appear. 

Two natural ways of simplifying Eq. (|I0|) are first integrating over uj and second setting 
uj = 0. These yield 

dcu ,,, . ij, B 



J -^(|/ l (q,^)| 2 ) = -{l + A 2 [ul 1 (d,z)-2ul 2 (d,z)]}\q\~ 

(\h(q, 0)| 2 ) = ^{1 + A + A 2 [1 + uMd, z) - AuJ 2 (d, *)]}|q| 
respectively, where r = d + 4 — 2^, u = 2g 2 B/A 3 , 

h(d,z) = J 
I 2 (d,z) = 



(12) 



(13) 



d d p 
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(P + l) 


1 2 

IP 


" r |P + l 


-\- T 


(27l) d 


fl + IPl 


z + 


P+l| z ) 





d d p P 1 p ■ (p + 1) 



(14) 
(15) 



(2vr) d (l + |p^ + |pj + i|*) ' 

and Ji(d,z) is the same as Ii(d,z) (i=l,2), except that the term in the denominator (1 + 
IpT + Ip + l| z ) is replaced by (|p| z + |p + l| z ). 



When carrying out the PMS scheme we ask that the correlation functions (12]) and 
(p!3|) be insensitive to the artificially introduced quantities (Tll|). Since z/ q and D q were 
assumed to have power-law forms, they are, apart from the amplitudes A and B, completely 
characterized by the exponent z. Thus, implementing the PMS in this case entails requiring 
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that the expressions (|P2|) and (|13D be independent of z. There is, however, a subtlety here. 
One cannot apply the PMS to a function and vary with respect to a single parameter as 
that results in an overdetermination. Either we perform a functional variation, or we apply 
the method to a q-independent quantity. Here we choose to vary with respect to z alone, 
so we should be analyzing a correlation function which has a q-independent scaling form. 
One possibility is to treat correlation functions of the type ((V a h) 2 ), which in q-space reads 
|q| 2a (|/i(q, cij)| 2 ), rather than (\h(q, u))\ 2 ). Here we choose a so that the analogues of (|12"D 
and fll3f) are q-independent. Differentiating these expressions with respect to z gives 

Scheme 1 d z [h(d, z) - 2I 2 (d, z)\ =0 (16) 

Scheme 2 d z [Ji{d, z) - U 2 (d, z)} = 0, (17) 

respectively. 

An alternative to starting from the full correlation function (|/i(q, w)| 2 ), is to extract the 
parts of this function which contribute to the renormalization of the noise strength. This is 
achieved by amputating the two external legs on the graphs for the correlation function, that 



is, removing two factors of G(q, u;) from ([TO). The resulting forms for the leading scaling 



behavior resemble (|I2|) and flT5|), but without the terms involving the integrals I 2 and J 2 . 
Applying the PMS to the noise strength gives equations analogous to fllBD and fll7|) : 

Scheme 3 d z [h(d, z)} = (18) 

Scheme 4 d z [Jx(d, z)\ = 0. (19) 

To find z from (p!6|)-([T9|) we have to evaluate integrals such as (|TJ]) and (|T5|) for a range 
of values of z. Going over to <i-dimensional spherical polar coordinates, the integrals may be 
reduced to double integrals depending on the parameters d and z. We have taken great care 
in numerically integrating the double integrals which give dli/dz and dJi/dz, (i = 1,2), 
since there are potential singularities at p = 0, — 1 and as |p| — > oo. In the first two cases 
the integrals converge if z > 1, which is always true. The singularities at large |p| are more 
severe: integrals of the type (0) only exist if z < (d + 4)/3 and those of the type ( |T5D if 



z < 2. In the latter case a naive analysis might suggest that the integral only exists when 
z < 1, but the leading large |p| term vanishes by p — > — p symmetry. 

Based on these considerations, we have searched for solutions of (pl^-(|T9|) for parameter 
values 1 < z < (d + 4)/3 for 1 < d < 2 and 1 < z < 2 for d > 2. Within this parameter 
range there are integrable singularities which numerical integration routines find difficult 
to cope with. For this reason, we have used various transformations which remove these 
integrable singularities, as well as different routines, in order to check our computations. 
Note that while schemes 3 and 4 seem to admit solutions with z > 2, it is easy to check 



that the subleading non-scaling terms which were discarded to obtain ( fL2|) and ( |13D are only 
subdominant if z < 2. Thus the entire method is dependent on this bound being respected. 

For Schemes 3 and 4, we readily find solutions, and the resulting values of z are reported 
in the table below. We find no solutions to Schemes 1 and 2 as presented above. However, 
the PMS methodology asserts that if the point of minimal sensitivity is not found by setting 
the derivative equal to zero, one should next seek to minimize the derivative, and so on ||. 
In this light we revise these schemes and look at higher derivatives. Scheme 2 has a solution 
when the second derivative is set to zero and the corresponding values of z are reported in 
the table. On the other hand, we have seen no evidence for a solution to Scheme 1 (in d — 2) 
upon examining the second and third derivatives. 

The simulation results included in the table are from the hypercube-stacking model 
simulated in dimensions d — 1,2,3 Jl3|] and the restricted solid-on-solid model, simulated 



up to dimension d = 7 [|T4|, both believed to belong to the KPZ universality class. The 
measurements of the exponent /3, which governs the early growth, agree between the models 
for d = 1,2,3. In the table we list the corresponding value for the dynamic exponent 
z = 2/(/3 + 1) and include the larger reported error bar. The data for d — 4 showed 
large fluctuations and finite-size dependence, making it more difficult to establish accurate 
estimates of the exponents ||14j| . 

Schwartz and Edwards adopted the FAC approach, determining z by setting the terms 
of 0(A) and higher in ([12]) to zero; the outcome of this procedure is included in the column 
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labeled SE. The column labeled BC contains the findings of Bouchaud and Cates who 
pursued a self-consistent approach. Actually the latter can also be formulated using the 
FAC criterion; the only difference being that SE integrates over u while BC sets uj = 0. 
Incidentally, the equations obtained by BC can also be derived as part of a functional 
PMS variation; but one in which scaling is still ultimately put in by hand, and not deduced. 

The values for z furnished by the PMS Schemes (2-4) follow the same general pattern: 
they fall below the exact value of z = 3/2 in d — lbya few percent; they then increase 
as d increases and eventually reach z = 2 in some critical dimension d c (listed in the table) 
in disagreement with the simulation results. Note that the self-consistent methods contain 
enough of the correct structure to yield the exact value of z = 3/2 in d = 1. The PMS 
approach uses this same structure and information but in a very different way and the exact 
result in d = 1 is no longer guaranteed. However, for d > 2, the PMS schemes are more in 
line with the simulation results. 

In this Letter we have applied the PMS to obtain improved estimates for the exponent 
z in the KPZ equation. The ideas discussed here could be extended in a number of ways. 
Besides the obvious, but tedious, course of going to higher order in the perturbation, one 
could adopt a more general (though still Gaussian) form for (^) allowing for a frequency- 
dependent noise strength or for a more general functional form G _1 (k, uS) in place of 
the response u^ — iu. This has proved successful within the self-consistent approach, where 
a more sophisticated ansatz for the response function in Ref. f| led to improved results 
for z. It would be interesting to see if more general ansatze such as these would similarly 
improve on the estimates presented here. A more ambitious calculation would be to apply 
the PMS to the function ( |10D without assuming scaling. In the absence of a renormalization 
group treatment of this problem, approaches such as these may be the most promising ways 
toward an understanding of the anomalous large- distance behavior of the KPZ equation. 
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Table 



d 


Simulation 


SE 


BC 


Scheme 2 


Scheme 3 


Scheme 4 


1 


1.500±0.001 


1.500 


1.500 


1.424 


1.420 


1.465 


2 


1.613±0.003 


1.705 


1.667 


1.582 


1.553 


1.621 


3 


1.695±0.007 


1.920 


1.862 


1.740 


1.682 


1.769 


4 


1.77±0.02 






1.909 


1.810 


1.914 


d c 




3.2 


3.6 


4.4 


5.4 


4.6 



Values of the dynamic exponent z. 
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